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Regular (non-singular) space-times are given which describe the formation of a (locally defined) 
black hole from an initial vacuum region, its quiescence as a static region, and its subsequent 
evaporation to a vacuum region. The static region is Bardeen-like, supported by finite density and 
pressures, vanishing rapidly at large radius and behaving as a cosmological constant at small radius. 
The dynamic regions are Vaidya-like, with ingoing radiation of positive energy flux during collapse 
and negative energy flux during evaporation, the latter balanced by outgoing radiation of positive 
energy flux and a surface pressure at a pair creation surface. The black hole consists of a compact 
space-time region of trapped surfaces, with inner and outer boundaries which join circularly as a 
single smooth trapping horizon. 

PACS numbers: 04.70.-s, 04.20.Dw 



Introduction. Black holes, predicted by Einstein grav- 
ity, appear to exist in the universe. The singularities 
which were also predicted to form inside them how- 
ever, are generally regarded as indicating the breakdown 
of the theory, requiring modifications which presumably 
include quantum theory. A first step in this direction, 
quantum field theory on a stationary black-hole back- 
ground, predicted Hawking radiation 2\. The ingoing 
radiation has negative energy flux which contradicts the 
assumptions of the singularity theorems and, in a semi- 
classical approximation, causes the black hole to shrink. 
In the usual picture [3|, the black hole shrinks until the 
central singularity is reached. However, if the singularity 
does not exist, such a picture cannot be correct 0,0. 

Regular (i.e. non-singular) black holes have sometimes 
been considered, dating back at least to Bardeen One 
can find metrics which are spherically symmetric, static, 
asymptotically flat, have regular centres, and for which 
the resulting Einstein tensor is physically reasonable, sat- 
isfying the weak energy condition and having components 
which are bounded and fall off appropriately at large dis- 
tance. The simplest causal structure is similar to that 
of a Reissner-Nordstrom black hole, with the internal 
singularities replaced by regular centres (Fig. ^ . Such 
space-times have been dismissed as unphysical, due to 
the presence of a Cauchy horizon, but if such a black 
hole evaporates, the Cauchy horizon is no more real than 
the event horizon, as examples will show. 

Imagine removing astrophysically irrelevant regions to 
the past and future of two consecutive advanced times 
(Fig. 2J, then adjoining a past which describes gravita- 
tional collapse and a future which describes evaporation. 
The static region contains inner and outer horizons which 
no longer have global significance, but still have local sig- 
nificance as trapping horizons 7J . The key issue is how 
the trapping horizons develop, which has been predicted 
on general principles 5j. In this Letter, concrete models 
are given for the collapse and evaporation phases, using 
Vaidya-like regions Q with ingoing or outgoing radiation. 

Regular static black holes. Consider static, spherically 




FIG. 1: Penrose diagram of a regular non-extreme black hole, 
to be identified vertically with isometric regions. 

symmetric metrics of the form 

ds 2 = r 2 dS 2 + dr 2 /F(r) - F(r)dt 2 (1) 

where t is the static time, r the area radius and dS 2 = 
dO 2 + dcj) 2 sin 2 8. A surface has area 47rr 2 , is trapped if 
F(r) < and untrapped if F(r) > 0. Trapping horizons, 
in this case also Killing horizons, are located at the zeros 
F(r) — 0, and there is a standard procedure to match 
regions across such horizons 9J. For an asymptotically 
flat space-time with total mass m, 

F(r) ~ 1 — 2m/r as r — ► oo. (2) 

Similarly, flatness at the centre requires 

F(r) ~ 1 - r 2 /l 2 as r -> (3) 
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FIG. 2: The metric function F — g rr , whose sign determines 
gravitational trapping, for fixed core radius I and different 
total masses m. 



where I is a convenient encoding of the central energy 
density 3/87rZ 2 , assumed positive. A sketch of F(r) in- 
dicating where it might dip below zero (Fig. |3J) shows 
that there will be a range of parameters for which there 
is no black hole, and that the simplest black-hole cases 
will generically have an inner and outer Killing horizon, 
the two cases separated by an extreme black hole with 
degenerate Killing horizon. 

It can be shown that, for a metric g of the form ifTjl. iffll. 
the Einstein tensor has the cosmological-constant form 
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Thus there is an effective cosmological constant at small 
distances, with Hubble length I. Such behaviour has been 
proposed previously by Sakharov [lfj as the equation of 
state of matter at high density, and by Markov and others 
based on an upper limit on density or curvature, to be 
ultimately justified by a quantum theory of gravity. Since 
I gives the approximate length scale below which such 
effects dominate, one might expect / to be the Planck 
length or of the same order, though larger length scales 
are not excluded. 

A minimal model. For definiteness, take a particularly 
simple metric satisfying the above conditions: 



F(r) = 1 
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(5) 



where (I, m) are constants which will be assumed posi- 
tive. This is similar to the Bardeen black hole, reduces 
to the Schwarzschild solution for I = and is flat for 
to — 0. Poisson & Israel [12| derived an equivalent form 
of g rr = F (without fixing g u ) based on a simple relation 
between vacuum energy density and curvature. 

Elementary analysis of the zeros of F(r) reveals a crit- 
ical mass to* = (3-\/3/ 4)Z and radius r* = \fil such that, 
for r > 0, F(r) has no zeros if to < to*, one double zero 
at r — r* if to = m* , and two simple zeros at r = r± if 
to > to* (Fig. [2J . These cases therefore describe, respec- 
tively, a regular space-time with the same causal struc- 
ture as flat space-time, a regular extreme black hole with 
degenerate Killing horizon, and a regular non-extreme 
black hole with both outer and inner Killing horizons, 
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FIG. 3: Horizon mass-radius relation: a pair of horizons ap- 
pears when mass m exceeds critical mass m„. 



located at r + « 2m and r_ « I for to to* (Fig. QJ. 
The horizon radii r± determine the mass (Fig. |2J 
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Note the existence of a mass gap: such black holes cannot 
form with mass m < to*. Also, the inner horizon has 
radius r_ > I which is very close to I for all but the 
smallest masses. In this sense, the black- hole core has a 
universal structure. 

If the Einstein equation G = 8nT is used to interpret 
components of the energy tensor T, these metrics are sup- 
ported by density — T t , radial pressure and transverse 
pressure Tg = T? given by 
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They fall off very rapidly, 0(r 6 ), at large distance. In 
terms of the energy E defined by 
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one finds the energy density ~T\ — (3l 2 /2-!r)(E/r 3 ) 2 1 
proportional to the square of the curvature E/r 3 . Pois- 
son & Israel 0] assumed such proportionality as a prop- 
erty of vacuum energy density; then the component 
dE/dr = —4irr 2 Tf of the Einstein equation implies g rr 
equivalent to (JSJl. 

Adding radiation. Next rewrite the static space-times 
in terms of advanced time 



v = t + 
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so that 
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Now allow the mass to depend on advanced time, m(v), 
defining F(r,v) by the same expression (0. Then the 
density — T", radial pressure and transverse pressure 
T[j have the same form Q-JSJl, but there is now an ad- 
ditional independent component, radially ingoing energy 
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FIG. 4: A mass profile m(v) in advanced time v. 

flux (or energy-momentum density) T£ given by 

2r 4 m' 
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(12) 



where to' = dm/dv. This describes pure radiation, recov- 
ering the Vaidya solutions for I — and at large radius. 
In the Vaidya solutions, the ingoing radiation creates a 
central singularity, but in these models, the centre re- 
mains regular, with the same central energy density given 
by ©■ It seems that the effective cosmological constant 
protects the core. 

The ingoing energy flux is positive if m is increasing 
and negative if m is decreasing. A key point is that trap- 
ping horizons still occur where the invariant g rr = F(r, v) 
vanishes ■ Then one can apply the previous analysis to 
locate the trapping horizons in (v, r) coordinates param- 
eterized by to, given by m(r±) in © and a mass profile 
m(v); qualitatively, by inspecting Figs. I^Jand^J 

Ingoing radiation. One can now model formation and 
evaporation of a static black-hole region. Introduce six 
consecutive advanced times v a < Vb < ■ . ■ < Vf and 
consider smooth profiles of m(v), meaning m'(v) at least 
continuous, such that (Fig. 0} 
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m{v) = 
m'(v) > 
m{v) = mo > m* 
m'(v) < 
m(v) = 0. 



Then 



3t>& € (v a ,v c ) : m(v b ) = 
3v e e {vd,Vf) : m(v e ) = m* 



(13) 
(14) 
(15) 
(16) 
(17) 



(18) 
(19) 



These transition times mark the appearance and dis- 
appearance of a pair of trapping horizons: for v < 
and v > v e , there is no trapping horizon, while for 
Vf, < v < v e , there are outer and inner trapping horizons, 
in the sense of the author's local classification 0. These 
horizons join smoothly at the transitions and therefore 
unite as a single smooth trapping horizon enclosing a 
compact region of trapped surfaces (Fig. for r < r ). 

Outgoing radiation. Thus far, only the ingoing Hawk- 
ing radiation has been modelled, since outgoing radiation 
does not enter the equation of motion of the trapping 
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FIG. 5: Penrose diagram of formation and evaporation of a 
regular black hole in the given models. 



horizon; in terms of retarded time u, T vv and T uv enter, 
but T uu does not 0- Outgoing Hawking radiation will 
now be modelled by adapting an idea of Hiscock [l3| : 
select a certain radius ro > 2mo outside the black hole, 
and adopt the above negative-energy radiation only in- 
side that radius, balanced by outgoing positive-energy 
radiation outside that radius, with the same mass pro- 
file (Fig. 03). This is an idealized model of pair creation 
of ingoing particles with negative energy and outgoing 
particles with positive energy, locally conserving energy. 
In more detail, consider an outgoing Vaidya-like region 
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(20) 



with F(r,u) as before (0, with m replaced by a mass 
function n(u). Fix the zero point of the retarded time 
u so that r — ro corresponds to u = v. Now take the 
above model only for v < Vd lfl"3")l - l(l"5|l . For v > Vd, keep 
the profiles (|16(l (|17|l for r < ro, but for r > ro, take an 
outgoing Vaidya-like region with 



Vu < Vd : n{u) = too 



(21) 
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Vw > Vd : n(u) = m{u). 



(22) 



Then there is a static region with total mass mo for v > 
Vd, u < Vd, and a flat region for v > Vf, u > Vf. Since 
the ingoing and outgoing radiation has no net energy but 
a net outward momentum, one might expect the pair 
creation surface r — tq to have a surface layer with no 
surface energy density but surface tension r < 0. This is 
confirmed using the Israel formalism [T3 |. yielding 
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at r = ro, Vd < v < vf. 

The whole picture is given in Fig. Action begins 
at v = v a , a black hole begins to form at v = Vb, has 
collapsed completely at v = v c to a static state with mass 
mo, begins to deflate at v = Vd and eventually evaporates 
at w = w e , leaving flat space finally after v = Vf, u = Vf. 
There is no singularity and no event horizon. 

Remarks. A trapping horizon with both inner and 
outer sections typically develops in numerical simula- 
tions of binary black-hole coalescence, in analytical ex- 
amples of gravitational collapse such as Oppenheimer- 
Snyder collapse and according to general arguments 0] . 
A key point here is that the inner horizon never reaches 
the centre, where a singularity would form This is 
compatible with the classical singularity theorems 0, 
which make assumptions that are already not satisfied 
by a Bardeen black hole, such as the strong energy con- 
dition. The negative-energy nature of ingoing Hawking 
radiation shows that such theorems do not apply to a 
black hole that might someday begin to evaporate. 



In contrast to the usual picture [3j , the endpoint v = v e 
of evaporation, defined locally by the disappearance of 
trapped surfaces, occurs when the outer and inner sec- 
tions of the trapping horizon reunite. The subsequent 
timescale until the effective cessation of particle produc- 
tion at v = Vf can be expected to be of the same order as 
I. Another logical possibility is that the inner and outer 
horizons approach each other asymptotically, forming the 
horizon of an extreme black hole with m = m* , but such a 
delicately balanced situation would require justification. 

The possibility of a circular trapping horizon has, in 
fact, been conjectured before |16|. Since there is no event 
horizon, long accepted as the defining property of a black 
hole, it seems necessary to stress that the static region 
looks just like a black hole over timescales that can be 
arbitrarily long. Thus it should be regarded as a black 
hole by any practical definition, as in the local, dynamical 
paradigm for black holes in terms of trapping horizons 0. 
The non-existence of an event horizon for a black hole 
which eventually evaporates seems to have been recently 
accepted by its most influential proponent [l7j . 

Most discussions of black-hole evaporation mention a 
certain I-word, as a paradox, problem or puzzle. The 
above space-times, regular and with the causal structure 
of flat space-time, show that this word need not be men- 
tioned. To paraphrase an old gravitational adage: what 
goes in, must come out. 
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